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Abstract. In non-hadronic axion models, which have a tree-level axion-electron interaction, 
the Sun produces a strong axion flux by bremsstrahlung, Compton scattering, and axio- 
recombination, the "BCA processes." Based on a new calculation of this flux, including for 
the first time axio-recombination, we derive limits on the axion-electron Yukawa coupling 
Qae and axion-photon interaction strength Qa^y using the CAST phase-I data (vacuum phase). 
For rua < 10 meV/c^ we find ga-y Oae < 8.1 x 10~23 QeV"^ at 95% CL. We stress that a 
next-generation axion helioscope such as the proposed lAXO could push this sensitivity into 
a range beyond stellar energy-loss limits and test the hypothesis that white-dwarf cooling is 
dominated by axion emission. 
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1 Introduction and main results 

The CERN Axion Solar Telescope (CAST) [1] is a helioscope experiment aiming at the 
detection of axions and axion-like particles (ALPs) emitted from the Sun. The detection 
principle is based on the axion^ coupling to two photons, which triggers their conversion into 
photons of the same energy as they propagate through a transverse magnetic field [2, 3]. 
CAST has tracked the Sun in three different campaigns (2003-04 [4, 5], 2005-2006 [6] and 
2008 [7]) with a 9.26 m long, 9 Tesla strong, decommissioned LHC dipole test magnet while 
measuring the fiux of X-rays at the exits of both bores with four different low-background 
detectors. No excess counts were observed over the expected backgrounds, thus constraining 
the axion parameters, notably mass and couplings to photons [4~7] and nucleons [8, 9]. 
Such constraints require knowledge of the fiux of axions emitted from the Sun which can be 
computed very precisely because the solar interior is a tractable weakly-coupled plasma. 

The physics case of CAST was mainly focused on hadronic axions [10, 11] which were 
appealing as hot dark matter candidates [12]. Hadronic axion models are minimal in that 
the generic axion interactions with hadrons and photons derive from mixing with the pseu- 
doscalar mesons vr'', rj and rj'. The interactions with leptons arise at loop level [13] and are 
usually irrelevant. In hadronic models, the bulk of the solar axion fiux comes from Primakoff 
production + Q ^ a + Q [14-16], where Q is any charged particle. 

Recently, non-minimal axion models are receiving increasing attention [17-26]. From the 
theoretical point of view, axions are nowadays known to arise naturally in many extensions 
of the standard model that pursue some sort of unification, such as Grand Unified Theories 
or string theory. Indeed, the original hadronic KSVZ axion [10, 11] can be regarded as an 
exemplary toy model that contains only the essential ingredients to solve the strong CP 
problem [27], but often axions arising in completions of the standard model are not minimal 
in this sense. 

^Unless otherwise noted, the term "axion" henceforth includes both QCD axions and more general ALPs. 
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Figure 1. Solar axion flux on Earth for a typical DFSZ model with interaction strength to photons 
g^^ = 10-12 GeV-i and electrons gae = 10"", corresponding to fa = 0.85 x 10^ GeV [30]. The blue 
line corresponds to the PrimakofF flux and the red lines show the different components of the BCA 
flux: FF = free-free (bremsstrahlung), FB = free-bound (axio-recombination), and BE — bound- 
bound (axio-dcexcitation). The black line is the total flux. 

Non-hadronic axion models such as that of DFSZ [28, 29] have different and very in- 
teresting phenomenological consequences. Notably, they couple to electrons at tree level, 
and this opens axion-production channels in stars which are much more effective than the 
Primakoff process: electron-ion bremsstrahlung [e + I^e + I + a) [15, 31], electron-electron 
bremsstrahlung {e + e — )• e + e + a) [15], Compton (7-l-e — )• e + a) [32, 33], axio-recombination 
(e -|- / — 7- /~ -|- a) [30, 34-36] and, to a lesser extent, axio-deexcitation of ions (/* — )• I + a). 
Henceforth we shall refer to this set of reactions as BCA for its most relevant contributions 
from bremsstrahlung, Compton, and axio-recombination. Indeed, axions with gae ~ 10"^"^ 
might explain the longstanding anomaly in the cooling of white dwarfs (WD) [37], recently 
reinforced by updated studies of the period decrease of the pulsating white dwarfs G117- 
B15A [38, 39] and R548 [40] and the WD luminosity function [41-44]. One should note that 
this value is somewhat challenged by the constraint <^ac < 2.5 x 10~^^ imposed by the evolu- 
tion of red giant stars in globular clusters [45, 46]. These values of gae imply a DFSZ-axion 
decay constant fa ~ 10^ GeV (corresponding to an axion mass ma ~ meV). Such me V- mass 
axions have a wealth of other interesting phenomenological implications in the context of 
astrophysics, like the formation of a cosmic diffuse background of axions from core collapse 
supernova explosions [47] or neutron star cooling [48, 49]. In cosmology, the decay of relic 
axionic strings and domain walls produces a relevant cold dark matter population [50]. 

Still, besides new studies of red giant evolution — currently underway [51] — or features 
in massive star evolution such as [52] we know of no other way to test the WD cooling 
hypothesis that could rely solely on the axion-electron coupling. It appears thus that to 
assess the WD cooling hypothesis we must study processes involving other axion couplings 
as well. In this paper we will make use of the coupling to two photons. As already mentioned, 
the solar BCA flux of non-hadronic axions is generically much larger than that of hadronic 
models for the same value of fa and moreover, it has a different spectrum (see figure 1). 
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Figure 2 . CAST constraints on Qae x ga-y as a function of , assuming the solar emission is dominated 
by the BCA reactions wliicli involve only the electron coupling gae- 

Therefore, helioscopes are appealing to search for non-hadronic axions. Indeed, some of us 
have recently shown that a next generation helioscope [53] , such as the proposed International 
AXion Observatory (lAXO) [54] can test the WD cooling hypothesis down to very small 

COUphngS Qae ~ 10~^^. 

In the present paper we analyze the CAST data in search of non-hadronic axions and 
set new upper bounds on gae x Qa^y, the product of the electron coupling (responsible for the 
production in the Sun) and the two-photon coupling (responsible for the detection in CAST). 

Figure 2 shows our results when we assume that the Primakoff emission from the Sun is 
subdominant and therefore the solar flux is caused by the BCA reactions alone. Our analysis 
of CAST data then constrains 

5a7 X Qae < 8.1 X lO'^^ GeV"^ (95% CL) (1.1) 

at low masses < 10 meV — where the probability of axion-photon conversion in CAST 
becomes independent of the mass — and worsens as for higher masses. 

If we also include the Primakoff flux (which is unavoidable because it is produced by 
the same coupling involved in the CAST detection), the signal at CAST depends inde- 
pendently on three parameters: gae, and ttIq. However, for nia < 10 meV the detection is 
independent of mass and we can plot our results in the gae^Qa^y parameter space. In this low- 
mass range, phase-I of CAST gives the strongest constraints and thus we have focused only 
on this data set. Our analysis, based on a two-free parameter likelihood method is able to 
exclude the region above the thick black line in figure 3. For very small values of gae ^ 10~^^, 
the BCA flux is negligible and the CAST bound smoothly becomes ga^y < 0.88 x 10~^° GeV~^ 
as found in our previous study [5] where only Primakoff emission was assumed. However, for 
larger values of gae the BCA flux becomes dominant and we recover equation 1.1. 

Note that our bound relies on a simple calculation of the solar axion flux, for which we 
have taken a solar model unperturbed by axion emission. If gae or gay are very large, the large 
axion flux requires a modified internal structure of the Sun with larger nuclear reaction rates 
and higher temperature of the core. The most stringent constraint derives from the agreement 
between the predicted and observed solar boron neutrino flux [55], excluding the gray region 
depicted in flgure 3 (labeled Solar u). Thus our bound is completely self-consistent up to 
gae = 3 X 10~^^, in contrast to those solar axion searches utilising Bragg scattering [56-59] 
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Figure 3. Constraints on gae and ga-y for nia ^ 10 mcV. The region above the thick black line 
is excluded by CAST. The gray region is excluded by solar neutrino measurements. In the vertical 
orange band, axion emission strongly affects white dwarf cooling and the evolution of low-mass red 
giants; parameters to the right of this band are excluded. Likewise, helium-burning stars would be 
perceptibly affected in the horizontal blue band; parameters above it are excluded. 



which have more limited sensitivity [60], as well as to some other searches relying solely on 
Qae coupling [61-63]. 

In order to put our results into context we also show in figure 3 two color bands rep- 
resenting parameters where axion emission would have a strong impact on stellar evolution. 
In the vertical orange band of gae values, axion emission would strongly affect WD cool- 
ing [37-44] and delay helium ignition in low-mass red giants [45, 46]. The exact range of gae 
values that is consistently ruled in or ruled out by these arguments remains to be studied in 
detail, but for sure parameters to the right of this band are excluded. Within the horizontal 
blue band, axion Primakoff emission would strongly affect stars in the helium-burning phase. 
The upper edge of this band corresponds to the traditional horizontal-branch star limit, the 
remaining range represents a new argument concerning the blue-loop suppression during the 
helium-burning phase of massive stars [52]. 

The orange band cuts the CAST constraint in its horizontal part which corresponds to 
the Primakoff flux dominating the solar flux, but very close to the values gae ~ 10~^^ where 
the BCA flux starts to dominate. Therefore, CAST cannot shed any further light on the WD 
cooling hypothesis. However, a next-generation helioscope such as lAXO with its improved 
sensitivity to (7a^ will also benefit from the large BCA-emitted flux and will improve over the 
RG bound in part of the parameter space. In principle, the WD cooling hypothesis is then 
testable in a laboratory experiment. 

After having presented our results and main messages, the rest of the paper is devoted 
to elaborate on our definitions, assumptions, and analysis method. In section 2 we give a 
brief account of axion theory, we further examine the implications of our findings, and finally 
describe the solar axion flux, and in section 4 we present our new analysis after a summary 
of the experimental setup of CAST phase-I and its results. 
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2 Properties of axions and cixion-like particles 



For the purpose of the present paper we can parametrize an axion model with the lagrangian 
density 

>C = \{d,a){d^a) - \mla^ - ^-fF.^F^^a - g,,^^^^,^'^^, , (2.1) 



where a is the axion field, its mass, Ffj^i, and F^'^ are the electromagnetic field tensor 
and its dual, mg the electron mass, and ipe the electron field. The coupling constant ga'y 
has units of energy"^ while gae is a dimensionless Yukawa coupling. Particles featuring this 
type of lagrangians are often called axion-like particles (ALPs) and they appear as pseudo- 
Nambu-Goldstone bosons (pNGB), associated with a global shift symmetry a — ?• a + const, 
which is spontaneously broken at some high energy scale fa (sometimes called ALP decay 
constant) or stringy axions where fa corresponds to the string scale Mg. The shift symmetry 
is explicitly broken by some perturbing dynamics responsible for the mass term. 

In 1977 Peccei and Quinn proposed one such symmetry to solve the strong CP problem 
[64, 65] with the additional condition that it should be color anomalous. The resulting pNGB 
was called the axion [66, 67]. The axion receives its mass from chiral symmetry breaking 
after mixing with the pseudoscalar mesons through the color anomaly, its magnitude being 

m^f^ lO^GeV 
rua = — 7 — - 6 meV , (2.2) 

^ ~T~ ^ Ja J a 

where is the neutral pion mass and is the pion decay constant. For the ratio z = mu/m^ 
of up to down quark masses we use the canonical value z ~ 0.56 although the allowed range 
is ^; = 0.35-0.60 [68], however this only leads to minor uncertainties in our context. 

The axion has a model-independent contribution to its two-photon coupling coming 
from the above-mentioned mixing with mesons and can also have a model-dependent part if 
the PQ symmetry has the electromagnetic anomaly. The two contributions sum to 

where a is the fine-structure constant and E/N the ratio of the electromagnetic and color 
anomalies of the PQ symmetry. 

The coupling to electrons has a model-dependent contribution proportional to an 0(1) 
coefficient X(, arising only in non-hadronic axion models and a very small model-independent 
one induced at one-loop via the photon coupling, 

me 3a^ rrie f E ^ fa , A \ , . 

fa 47r fa \N me me J 

where A is an energy scale close to the QCD confinement scale. 

For a generic ALP, (p, we expect similar equations for g^e and g^j,^ as for the axion, of 
course after removing the terms coming from axion-meson mixing, (the terms involving z) 
and changing fa for the ALP decay constant, However, the ALP mass is then completely 
unrelated to the couplings. 
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3 Expected counting rate 



3.1 Solar flux 

Based on the different axion interactions, different processes contribute to the production of 
the solar flux of these particles (see figure 4). The most important processes are: 

• Primakoff effect: 7 + Q — )• Q + a 

• Compton scattering (photo production): 7 + e — )• e + a 

• Electron-Ion bremsstrahlung (free- free transition) : e + /— )-e + / + a 

• Electron-electron bremsstrahlung: e + e— T-e + e + a 

• Axio-recombination (free-bound transition): e + / — t- /~ + a 

• Axio-deexcitation (bound-bound transition): /* — t- / + o 

where Q stands for any charged particle in the solar plasma, e for electrons, / for ions and 
/* for their excited states. The Primakoff process depends on the two-photon coupling and 
dominates when the coupling to electrons is absent at tree-level. When this is not the case 
the BCA processes dominate: bremsstrahlung on hydrogen and helium nuclei dominates the 
emission of low-energy axions, axio-recombination of metals (mostly O, Ne, Si, S and Fe) 
contributes sizably at intermediate energies and Compton takes over at higher energies. The 
contribution of axio-deexcitation is dominated by Lyman transitions (mostly Ly— a) and 
is significant only in the case of iron which dominates the axion flux around ~ 6.5 keV. 
Following Ref. [15] and integrating the emission rates over a solar model [69], the following 
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Figure 4. Some Fcynman diagrams for the most relevant solar axion emission reactions included in 
this work. 

fits for the axion fluxes at Earth can be obtained [53] (in units of year"""^ keV~"^) 
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where u) is the axion energy in keV, while P, C and B stand for Primakoff, Compton, and 
bremsstrahlung, respectively. In particular, the bremsstrahlung flux includes both electron- 
nucleus (only H and He) and electron-electron contributions. As a novelty, we include the 
emission of axions in the electron capture by an ion, dubbed "axio-recombination" (free- 
bound transition) and atomic "axio-deexcitation" (bound-bound transition). Originally, the 
former process was estimated to be subdominant in the Sun [34, 35]. However, a new cal- 
culation to be presented elsewhere [30], including the missing factor of 2 in the cross-section 
pointed out in [36] and captures in higher shells than the K-shell, shows that these processes 
are significant and increase the total flux by a factor of order 1 (figure 1). Unfortunately, 
the kinematic edges and the narrow lines from bound-bound transitions prevent us from 
providing a simple fitting formula. 

3.2 Helioscope event number 

The expected number of photons N-y from axion conversion in a given detector is obtained 
by integrating the product of the differential axion fiux, conversion probability and detection 
efficiency over the total range of energies 

= / da;(— ^) Va^^eSt (3.4) 

Juio V duj / total 

where S is the detection area perpendicular to the flux of axions, t is the exposure time, and 
e the detection efficiency. The axion-photon conversion probability in a transverse homoge- 
neous magnetic field B over distance L is 

ft.-, = (^)'w(f). (3^5) 

where sincx = (sinx)/x and the momentum transfer provided by the magnetic field is g' = 
m}j2u3. Coherent a-7 conversion along the full magnetic length happens when < ^bjjL^ 
i.e. when the momentum transfer is smaller than about where sinc(gL/2) — )■ 1. For 

larger masses, the conversion is not coherent^ and the probability gets suppressed by a factor 
~ (4a;/m^L)^. This factor is responsible for the degradation of our bound for vtia > 10 meV 
seen in figure 2. 

4 CAST experiment and analysis 

The most sensitive helioscope to date is the CERN Axion Solar Telescope (CAST), which 
makes use of a prototype superconducting LHC dipole magnet providing a magnetic field of 
up to 9T. CAST is able to follow the Sun twice a day during sunrise and sunset for a total 
of 3h per day. At both ends of the 9.26 m long magnet. X-ray detectors [70-72] have been 
mounted to search for photons from Primakoff conversion. 

CAST began its operation in 2003 and, after two years of data taking, determined an 
upper hmit on < 0.88 x IQ-^o GeV"^ at 95 % CL for < 0.02 eV [4, 5] in the hadronic 
axion scenario. To extend the experimental sensitivity to larger axion masses, the conversion 
region of CAST was filled with a suitable buffer gas [73] providing the photons with an 
effective mass, yielding upper limits on 5^7 ^ 2.3x 10~^'^GeV~^ for 0.02 < ttIq < 0.65 eV [6, 7]. 
Higher axion masses have also been studied and future publications will report on these 
masses which cover a range up to 1.18 eV. 

^We only use data from the CAST vacuum phase where photon refraction by the residual gas can be 
neglected due to the high vacuum conditions [5]. 
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4.1 The X-ray telescope of CAST 

The most sensitive detector system operative at CAST phase-I was the X-ray telescope [72], 
a combination of X-ray mirror optics [74] and a Charge-Coupled Device (CCD) [75] located 
in the focal plane of the mirror optics. Both instruments were originally built for satellite 
space missions, and together they increase the axion discovery potential significantly along 
with providing excellent imaging capability. The implementation of the X-ray mirror optics 
suppresses background by a factor of ~ 155, since photons are focused from the magnet 
aperture area of about 14.5 cm^ to a spot of roughly 9.3 mm^ on the CCD chip. One of the 
resulting advantages of focusing optics is the possibility to measure background and signal 
simultaneously. 

4.2 Data taking 

The CCD detector showed a stable performance over the entire 2004 running period. For 
our analysis, we have used a total of 197 h of tracking data (i.e. magnet pointing to the 
Sun) and 1890 h of background data taken from the same area during non-tracking periods. 
For a detailed description of the X-ray telescope design, its performance, and background 
systematics during the 2004 data taking period we refer to Ref. [70]. Since no significant 
signal over background was detected with the X-ray telescope during the data taking period 
of 2004, upper limits for the non-hadronic axions were derived for this CAST detector system. 

7 
6 
5 

3 

o 

" 3 
2 
1 



" 1 2 3 4 5 6 

a [keV] 

Figure 5. CCD spectra for tracking (blue) and background (red) runs during 2004. In both cases, the 
error bars represent the statistical uncertainty of the measurement. Please notice that the background 
spectrum has been renormalized to the tracking time of 197 hours. 

In order to minimize the influence of the Cu-Kq fluorescence line (at 8 keV originating 
from the cooling mask of the detector [70]), we restricted our analysis to the energy range 
between 0.8 and 6.8 keV. In total, we observe 26 counts in this energy range inside the signal- 
spot area during axion sensitive conditions. The background, defined by the data taken from 
the same CCD spot region during non-tracking periods, has been acquired under the same 
operating conditions. The spectral distribution of the observed events during tracking and 
non-tracking times with the CCD detector is shown in figure 5. 
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4.3 Spectral fitting results 

The resulting low counting statistics required the use of a maximum-likelihood method to 
determine an upper limit on gaeda-y- T'^g likelihood function used is based on a Poissonian 
p.d.f., the binned likelihood 

3 ■' 

where n = 20 is the number of spectral bins, tj the number of observed counts in tracking, 
and Aj the value of the mean in the j-th bin, respectively. The fit function, Aj = aj + bj 
is used, where bj is the measured background and aj oc gaeQa'y the expected number of 
counts in the j^^ energy bin from axion-to-photon conversion. The best estimate for gaeda-y 
is obtained by minimizing = — 21n£. The validity of the x^-interpretation in our case, as 
well as the negligible influence of the statistical uncertainty of the background on the final 
result have been verified with a Monte Carlo model by means of the generation of pseudo-data 
sets. 

24.0 
23.5 
23.0 
22.5 
22.0 
21.5 
21.0 
20.5 
20.0 

'^■^ -0.006 -0.004 -0.002 0.000 0.002 0.004 0.006 

g2_^gyi0-«GeV-] 

Figure 6. ^s, function oi g^^g^^^ for an axion mass of 1 meV. The mimimum of the distribution, 
— 1.136 X lO^'^^ GeV^^, is the most probable value of gl^gl^- 

We compared the result derived with the maximum- likelihood defined in equation 4.1 
with a maximum-likelihood technique based on an unbinned maximum-likelihood estimator 
that divides the exposure time in sufficiently small time fragments so that either one or zero 
counts are found in the detector. The unbinned likelihood can be expressed as 

N 

logC^-RT + ^logRitk,ujk) (4.2) 

k 

where the sum runs over each of the N detected counts and R{tk,ujk) is the event rate at the 
time tfc, energy ojk of the /c-event. Rt is the expected number of counts over all exposure 
time and energy 

R{t,u:) = B + S{t,uj) (4.3) 
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where B is the background rate of the detector and S{t,uj) is the expected rate from axions 
(A/"^), which depends on the axion properties ga^y and rua (see equation 3.4). 

The final analysis yields an upper limit on Qaeda-y ^ 8.1 x 10~^^GeV~^ (95% CL) for 
axion masses nia < lOmeV (see figure 2). Both likelihood methods give compatible results 
on both, the best fit value and the upper limit. In the latter case the deviation of the limit 
remains within ~ 0.3% over the considered axion mass range. 

Parallel to the methods described above, a two-free parameter likelihood was applied for 
axion masses within the reach of CAST. In this case, not only is the axion-electron coupling 
taken into account but also the Primakoff contribution to the axion production in the Sun in 
the frame of non-hadronic models. This approach computes a three-dimensional probability 
function at a given axion mass that correlates both ^^-y and gae couplings. The result yields 
a limit on the ga^-gae parameter space (see figure 3). 

4.4 Systematic uncertainties 

We studied the infiuence of systematic uncertainties on the best fit value of gae9aj ^-nd on 
the upper limit. Statistically significant variations of the background on long and short time 
scales are not apparent for the X-ray telescope data [70]. Even if the background level were 
time-dependent, it would not play a significant role, since the X-ray telescope is measuring 
both, potential signal and background, simultaneously. Since we observed no significant 
difference between the background level and the spatial distribution during tracking and non- 
tracking times, we used the same signal-spot area during non-tracking periods to define the 
background. Alternatively, we used tracking data and different regions on the CCD outside 
the signal spot area to estimate the systematic uncertainties due to the choice of background 
definition. The overall systematic uncertainty is dominated by both the background definition 
and the pointing accuracy, that affects the effective area of the telescope, the location of the 
signal spot and its size, respectively. Other effects such as uncertainties in the detector 
calibration, magnet parameters and the likelihood method used are negligible in comparison 
with the systematic induced by the choice of background definition. For the best fit value of 
9ae9a-y '^^ find that in the axion mass range for which CAST remains coherent we obtain 

glal^lhcstm = (-1.136 ±lil Stat. ±i|^ syst.) x 10^^^ GeV-^ . (4.4) 

Alternatively, the background was also determined by extrapolating the background 
measured during tracking periods in the part of the CCD not containing the Sun spot. 
Different background selections led to different upper limits on gaeda^/^ of them within the 
statistical uncertainty. 

5 Conclusions 

Axions with tree-level coupling to electrons provide a different physics case and phenomenol- 
ogy than hadronic models. From the theoretical point of view, non-hadronic models are 
appealing because they arise in grand unified theories (GUTs), well-motivated completions 
of the standard model at high energies. From the phenomenological side, it is worth noting 
that naturally the coupling to electrons leads to larger axion fiuxes from stars than the cou- 
pling to photons. The sensitivity of CAST to non-hadronic axions allows us to set a bound 
on the product of both coupling constants gaeOa^y ^ 8-1 x 10"^'^ GeV~^ for nia < lOmeV. 
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For hadronic axions, Primakoff emission is the dominant production process in stars 
and helioscope hmits depend only on the axion-photon interaction strength ga^. For low- 
mass axions, the CAST limit on ga'y is competitive with, and even somewhat superior to, 
the energy- loss limits from globular cluster stars. For non-hadronic models, the stellar fluxes 
are dominated by the BCA processes that are based on the axion-electron coupling Qae- In 
this case CAST constrains the product gaeda'y in a significant way. However, the stellar 
energy-loss limits on gae are so restrictive here that CAST is not yet competitive. 

The claim of an anomalous energy loss in white dwarfs could be an indication for the 
existence of axions with a coupling to electrons in the gae ~ 10~^^ range, close to the energy- 
loss limits of red giants in globular clusters. If we are to assume this hint, the flux of solar 
axions is fixed by this parameter and a next generation axion helioscope such as lAXO could 
be able to detect it. We believe that the test of the white-dwarf claim and surveying the DFSZ 
parameter space for the first time in a laboratory experiment is a compelling motivation for 
lAXO and strengthens its physics case. 

In summary, there is a strong motivation to improve the helioscope sensitivity beyond 
CAST down to ga'y ~ 10~^^ GeV~^ and gae ~ 10~^^ [53]. This region includes, on the 
high-mass end, a large set of favored QCD axion models, potentially supersedes limits from 
SN 1987A and red giants in non-hadronic models, and starts probing the parameters sug- 
gested by white-dwarf cooling. 
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